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Summary Background: Human coronavirus (HCoV) OC43 is the most prevalent HCoV in respi¬ 
ratory tract infections. Its molecular epidemiological characterization, particularly the geno- 
typing, was poorly addressed. 

Methods: The full-length spike (S), RNA-dependent RNA polymerase (RdRp), and nucleocapsid 
(N) genes were amplified from each respiratory sample collected from 65 HCoV-OC43-positive 
patients between 2005 and 2012. Genotypes were determined by phylogenetic analysis. 
Recombination was analyzed based on full-length viral genome sequences. Clinical manifesta¬ 
tions of each HCoV genotype infection were compared by reviewing clinical records. 

Results: Sixty of these 65 samples belong to genotypes B, C and D. The remaining five strains had 
incongruent positions in the phylogenetic trees of the S, RdRp and N genes, suggesting a novel 
genotype emerging, designated as genotype E. Whole genome sequencing and bootscan analysis 
indicated that genotype E is generated by recombination between genotypes B, C and D. Tem¬ 
poral analysis revealed a sequential genotype replacement of C, B, D and E over the study period 
with genotype D being the dominant genotype since 2007. The novel genotype E was only de¬ 
tected in children younger than three years suffering from lower respiratory tract infections. 
Conclusions: Our results suggest that HCoV-OC43 genotypes are evolving. Such genotype shift 
may be an adapting mechanism for HCoV-OC43 maintaining its epidemic. 
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Introduction 

Coronaviruses (CoVs), belonging to the family Coronaviri- 
nae, are a large group of viruses with a broad infection 
spectrum in human and animals. CoVs are related to respi¬ 
ratory tract disorders, gastroenteritis, as well as to sys¬ 
temic and neurological diseases. CoVs are the largest 
RNA viruses, containing a positive-sense, single-stranded 
RNA genome with a length of 27,000—31,500 nucletides. 
Based on genome phylogeny and serological characteriza¬ 
tion, CoVs are divided into four genera, Alphacoronavirus 
(a-CoV), Betacoronavirus (p-CoV), Gammacoronavirus (y- 
CoV), and Deltacoror\avirus(h-CcN). Since the isolation 
of HCoV-229E and -OC43 in 1960s, a total of six HCoV spe¬ 
cies have been identified, including severe acute respira¬ 
tory syndrome CoV (SARS-CoV) in 2003, NL63 and HKU1 in 
2004, and middle east respiratory syndrome CoV (MERS- 
CoV) in 2012. 1,4 HCoVs belong to a- (229E and NL63) and 
p-genera (OC43, HKU1, SARS-CoV and MERS-CoV). 

HCoVs were previously not considered to be of great 
importance with respect to human diseases as most HCoV- 
infections were thought to be associated with mild symp¬ 
toms and occasional lower respiratory tract infections 
(LRTIs) until an outbreak of SARS in 2003. That has led to 
increased concerns about HCoVs, while the identification of 
MERS-CoV in 2012 reinforced the public health significance 
of HCoVs. Although SARS-CoV is no longer detected since 
2004, MERS-CoV continued as an epidemic, spreading to 
more patients and countries. This spread indicates a high 
adaption capability of MERS-CoV in humans. Insight into 
the epidemic characteristics of HCoVs at the molecular 
level will allow us to predict viral pathogenesis and trans¬ 
mission activities and inform HCoV prevention and control, 
particularly against newly emerging HCoVs. 

HCoV-OC43 has been more prevalent than other common 
HCoVs including HCoV-229E, —NL63 and —HKU1, in pediat¬ 
ric and adult respiratory infections, and can also cause 
outbreaks in human respiratory tract infections. How¬ 
ever, our understanding of the molecular epidemiology of 
HCoV-OC43 has been very limited. The genetic diversity 
of HCoV-OC43 was first reported in Belgium in 2005 and 
three clusters were identified based on the analysis of the 
spike (S) gene of the prototype strain ATCC VR-759 and 
seven clinical strains. Subsequently, Lau et al. gave the 
first description on the molecular epidemiology of HCoV- 
OC43 using sequences from 29 clinical samples in 2011. 
Four genotypes, A, B, C and D, were identified based on 
the viral genome and the phylogeny of the main structural 
genes, S, RNA-dependent RNA polymerase (RdRp), and 
nucleocapsid (N) genes, and genotype D was reported to 
have arisen due to natural recombination. 1 However, these 
observations were based on only a limited number of HCoV- 
OC43 positive cases. Due to the limited availability of virus 
sequences, the molecular epidemiological characterization 
of HCoV-OC43, particularly its genotyping, was poorly 
deciphered. 

In this study, we genotyped HCoV-OC43 by analyzing full- 
length sequences of S, RdRp, N genes and viral genomes 
directly from respiratory samples collected from 65 HCoV- 
OC43 positive patients with acute respiratory tract in¬ 
fections (ARTIs) recruited from 2005 to 2012. We observed 


a genotype shift in HCoV-OC43 over the study period and 
confirmed the emergence of a new genotype E arising 
through natural recombination. 

Methods 

Patients and clinical specimens 

Patients suffering from ARTIs were recruited from the 
Beijing Children Hospital and the Peking Union Medical 
College Hospital in Beijing, China from March 2005 to 
December 2012 when they seek health care at these 
hospitals. Criteria for including patients in our study 
encompassed acute fever (body temperature >37.5 °C) 
with respiratory symptoms such as cough or wheezing, 
normal or low leukocyte count, and with or without 
radiological pulmonary abnormalities. Nasopharyngeal as¬ 
pirates (NPAs) were collected from pediatric patients. 
Nasal and throat swabs were collected from adult patients. 
The respiratory samples were stored in viral transport 
medium (VTM) at -80 °C before use. Clinical information 
of each enrolled patient was recorded in standard form and 
reviewed retrospectively. Written informed consent was 
obtained from all participants or guardians on behalf of the 
minors/children participants. The study was approved by 
the Medical Ethic Review Board of the Institute of Pathogen 
Biology, Chinese Academy of Medical Sciences. 

Molecular detection of HCoVs 

Viral nucleic acids were extracted from 200 pi respiratory 
samples using a NucliSens easyMAG apparatus (bioMerieux, 
Marcy I’Etoile, France) according to the manufacturer’s 
instructions and were stored at -80 °C until use. HCoV- 
OC43 positive respiratory samples were tested by RT-PCR 
with HCoV-conserved primers and were confirmed by 
sequencing methods as described elsewhere. The pres¬ 
ence of other common respiratory viruses was also deter¬ 
mined as described elsewhere, including influenza virus 
(IFV) A, B and C, human parainfluenza virus (HPIV) 1—4, 
adenovirus (Adv), respiratory syncytial virus (RSV) A and 
B, human metapneumovirus (hMPV), human bocavirus 
(HBoV), rhinovirus (HRV) and enterovirus (HEV). 14 

Sequencing of HCoV-OC43 genes and viral genome 

Total RNA from respiratory specimens was converted to 
cDNA using combined random primers and oligo(dT) primers 
and the Superscript III reverse transcription system (Invi- 
trogen, Carlsbad, CA). The full-length S, RdRp, N genes and 
viral genomes were amplified from each respiratory spec¬ 
imen which was positive for HCoV-OC43, using specific 
primers (Table SI) with a genome walking method. PCR was 
performed using the following conditions: 94 °C for 5 min, 
40 cycles of amplification at 94 °C for 30 s, 50 °C for 30 s, 
and 72 °C for 90 s, with a terminal elongation step at 72 °C 
for 10 min. PCR products were sequenced directly using an 
ABI 3700 DNA sequencer (Applied Biosystems, USA). Se¬ 
quences were assembled manually through alignment to 
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the reference strain HK04-02 (GenBank accession no. 
JN129835). 

Phylogenetic analysis 

All HCoV-OC43 sequences available in GenBank (www.ncbi. 
nlm.nih.gov) were retrieved on May 30, 2013. The back¬ 
ground information of all the sequences used for phyloge¬ 
netic analysis is summarized in Table S2. The full-length 
S, RdRp, N genes, and viral genomes of HCoV-OC43 were 
aligned using ClustalW program implemented in MEGA 5.1 
with sequences deposited in GenBank. Pair-wise 
sequence identities in each region were calculated for 
the comparison of sequence divergence using BioEdit. 
Maximum likelihood (ML) trees were constructed with the 
best fit model of General Time Reversible with gamma- 
distributed rate variation across sites and 1000 bootstrap 
pseudo-replicates implemented in MEGA 5.1. The bovine 
coronavirus was used as the outgroup sequence, but is 
not shown in the presented figures to make the phyloge¬ 
netic relationships more clear. Substitution models were 
selected using Modeltest (version 3.7) according to the 
Akaike information criterion. Phylogenetic trees of each 
gene region of HCoV-OC43 were constructed by using the 
neighbor-joining method with Kimura’s two-parameter 
model and 1000 bootstrap pseudo-replicates implemented 
in MEGA 5.1. 3 To analyze the recombination events, the 
genomes of HCoV-OC43 were aligned and analyzed using 
boot scanning method implemented in SimPlot (V3.5.1, 
http://sray.med.som.jhmi.edu/SCRoftware). 

Statistical analysis 

Distribution frequencies of HCoV-OC43 genotypes were 
compared by using Pearson’s Chi square test or Fisher’s 
exact test. One-way analysis of variance was used to 
analyze the continuous variables for population parame¬ 
ters. P values <0.05 were considered statistically 
significant. 

Nucleotide sequence accession numbers 

The nucleotide sequence data of S, RdRp, N genes and viral 
genomes of HCoV-OC43 used in this study have been lodged 
in GenBank and the accession numbers are shown in 
Table S2. 

Results 

Genotyping of HCoV-OC43 strains 

To genotype the HCoV-OC43 samples, we constructed ML 
trees using the full-length sequences of S, RdRp and N genes 
amplified from the 65 respiratory samples of HCoV-OC43 
positive patients in this study and compared them to those 
retrieved from GenBank (Fig. 1). The HCoV-OC43 sequences 
fell into four distinct clusters on the phylogenetic tree of 
the S gene as reported by Lau et al. 1 However, incongru¬ 
ities were observed in ML trees of RdRp and N genes, indi¬ 
cating genetic diversity. Briefly, OC43 strains identified in 


this study (designated CN strains) fell into three clusters 
in S gene, i.e., B, C and D genotypes, similar to those 
from Hong Kong, China (HK) and Belgium (BE). Eleven CN 
strains fell into genotype B together with five 2004 HK 
strains and the Belgium strain BE03. The sequences of 
this genotype possessed nucleotide (nt) identities of 
98.7%—99.6%. Three CN strains and 15 HK strains formed ge¬ 
notype C, possessing 99.6%—99.8% nt identities; while 51 CN 
strains clustered with nine HK strains and a BE04 strain to 
form genotype D, possessing 99.3%—100% nt identities. Ge¬ 
notype A contained only the cell culture strain ATCC VR-759 
as previously reported. 

The strains that fell into genotype C clustered together 
in the ML tree of the RdRp gene, as well as in the ML tree of 
the S gene. Strains belonging to genotype D clustered 
together with strains of genotype B, and these sequences 
possessed 99.7%—99.8% nt identities. Notably, five CN 
strains (1783A/10, 2058A/10, 2941 A/11, 3074A/12 and 
3194A/12), which belong to genotype B in the ML tree of 
the S gene, formed a distinct clade in the tree of RdRp 
gene. Multiple alignment of RdRp results showed that these 
five CN strains possessed 99.5—99.6% nt identities to 
B_BE03, CJHK04-01 and DJHK11-01, while other B strains 
possessed 99.7—100% nt identities to B_BE03 (Table 1). 

Analysis of the N genes showed that the strains that 
belong to genotype B (other than the five distinct CN stains) 
in the ML tree of the S gene clustered together, while the 
strains belonged to genotype C and D in the ML tree of the S 
gene clustered together. The aforementioned five distinct 
CN strains were separated from all the known genotypes 
and formed two clades. Multiple alignment results were 
consistent with our phylogenetic analysis as the five distinct 
CN strains had lower nt identities with representatives of B, 
C and D genotypes than other genotype B strains had with 
the reference strain, including B_BE03 (97.6-98.7%), 
CJHK04-01 (97.6-99.1%) and DJHK11-01 (97.5-99.0%) 
(Table 1). Taken together, the incongruities in the phyloge¬ 
netic trees together with the analysis of nt identities 
showed that a novel genotype, may have arisen, which 
we designated as genotype E. 

Recombination analysis of genotype E strains 

The incongruent phylogenetic pattern of the S, RdRp and N 
genes in the five genotype E strains, particularly the drop¬ 
out of 1783A/10 from the linage formed by other genotype 
E strains in the phylogenetic tree of N genes, indicate the 
occurrence of potential recombination events. To further 
demonstrate the emergence of genotype E strains, we 
amplified the whole viral genome sequences directly from 
respiratory samples. We obtained the whole genome 
sequences of four of the five distinct strains (1783A/10, 
2058A/10, 3074A/12 and 3194A/12; 2194A/11 was not 
available due to the very low viral load in the specimen). 
We then analyzed the potential recombination by con¬ 
structing the phylogenetic trees of all known 23 gene 
regions of these four strains. Ten other whole genome 
sequences of OC43 were used as reference strains, 
including BE03 and 2145A/10 (genotype B), HK04-01 and 
3647/06 (genotype C), BE04, HK04-02 and 5240/07 (geno¬ 
type D), and the ATCC strain (genotype A) (Fig. 2). Bovine 
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Figure 1 Phylogenetic analysis of different HCoV-OC43 strains based on full-length S, RdRp and N genes. Trees were constructed 
using the maximum-likelihood method. Strains identified in this study are presented in bold. A, B, C, D and E represent genotypes. 


CoV (accession no. U00735) was used as outgroup sequence, 
which was not displayed in the figure to save spaces. We 
found that these four strains form a separate linage (geno¬ 
type E) in the phylogenetic trees of complete genome, S, 
RdRp and most of the nonstructual protein (nsp) genes. 
These findings further confirmed that these distinct CN 
strains belong to a novel genotype E, despite the incon- 
gruent phylogenetic pattern was observed in ns5a, E, M 
and N genes. 

Notably, in the phylogenetic trees of the nsp2-nsp6 
genes, these four genotype E strains were closely related 


to genotype C; while clustered more closely with the strains 
of genotype B in the trees of nspl, nsp8, hemagglutinin- 
esterase (HE) and the S genes. Strains 3074A/12 and 3194A/ 
12 were also clustered together with genotype D in enve¬ 
lope (E) and membrane (M) genes. These results support our 
hypothesis that recombination events occur among OC43 
genotypes. 

To verify these findings, we then carried out boot 
scanning analysis and the genome sequences of B_2145A/ 
10, C_3647/06 and D_5240/07 were used as references. 
When the genomes of 1783A/10, 2058A/10, 3074A/12 and 
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Table 1 Nucleotide sequence identities for the S, RdRp and N genes of genotype E compared with the reference sequences of 
genotypes B, C and D. 


Phylogenetic trees by gene 

Genotype 

Reference strains of each genotype 

B_BE03 C_HK04-01 

D_HK11 -01 

S 

E 

98.7—99.1 a 

96.6-97.2 

96.0-96.7 


B b 

98.8-100 

96.6-97.2 

96.0-96.7 

RdRp 

E 

99.6 

99.5-99.6 

99.5-99.6 


B 

99.7-100 

99.4-99.6 

99.7-99.8 

N 

E 

97.6-98.7 

97.6-99.1 

97.5-99.0 


B 

99.4-100 

99.1-99.5 

99.1-99.4 


BE03: 87309 Belgium 2003, GenBank accession no. AY903459. 
HK04-01: GenBank accession no. JN129835. 

HK11-01: GenBank accession no. not available. 
a Nucleotide identity. 

b Sequences belong to genotype B in S, RdRp and N genes. 


3194A/12 were used as query sequences, we identified 
several potential recombination sites in the viral genomes 
of genotype E (Fig. 3). Here 3074A/12 was used as an 
example to show the recombination analysis results. From 
positions nt 1000 to 14,500, most of the region of 3074A/ 
12 were closely related to C_3647/06, except positions up¬ 
stream of nt 1,000, nt 2500 to 4,500, and nt 11,500 to 
12,500, where 3074A/12 was closely related to B_2145A/ 
10. From positions of nt 14,500 to nt 28,000, most of the re¬ 
gion was closely related to B_2145A/10. From positions nt 
28,000 to the 3' end of the viral genome, most of the region 
was closely related to D_5240/07. Potential recombination 
sites were at the junctions of nsp2/nsp3, nsp6/nsp7, nsp9/ 
nspIO, nsp12/nsp13, ns5a/E and M/N corresponding to the 
schematic diagram of the whole viral genome (Fig. 3). 
These findings were consistent with the observations in 
phylogenetic analysis of S, RdRp and N genes described 
above. Similar boot-scanning results were obtained when 
3194/12 was used as query strain. Most of the recombina¬ 
tion sites were also found when 1783A/10 and 2058A/10 
were used as query strains. However, lower similarities 
were found in ns5a, M and N gene regions between 
1783A/10, 2058A/10 sequences and references, which indi¬ 
cates the diversity of parent strains of recombination. 

Taken together, these findings indicate that natural 
recombination events led to the emergence of novel 
genotype E and suggest complicated recombination events 
in the circulation of HCoV-OC43 strains in nature. 

Temporal evolution of HCoV-OC43 genotypes 

Genotype shift plays an important role in virus adaption to 
hosts. To determine whether genotype shift occurred 
in HCoV-OC43, the yearly distribution of genotypes during 
the study period (2005—2012) were determined. HCoV- 
OC43 positive cases were identified for each year analyzed, 
and their detection rate ranged from 1.9%„ to 13.9% 0 with 
the highest detection rates in 2007 (Fig. 4). We found 
that the detection rate of HCoV-OC43 spiked every other 
year except in 2010. Shifts of HCoV-OC43 genotypes over 
time were observed. After a low level epidemic of geno¬ 
types C and B, genotype D became the major epidemic 
since 2007, with the highest detection rate of HCoV-OC43 


during the study period, and dominated between 2007 
and 2009. Genotype B re-occurred in 2010 together with 
the novel genotype E. The co-epidemics of genotype B, D 
and E were observed in 2010 and 2012. Genotype C has 
not been detected since 2006. 

Clinical characteristics of HCoV-OC43 genotype 
infections 

To characterize the clinical manifestations of different 
HCoV-OC43 genotypes, the clinical data of the 65 HCoV- 
OC43 positive cases were analyzed (Detailed information of 
each patient is summarized in Table S3). Of the 65 cases, 28 
were children less than 14 years old, one was a 16-year-old 
teenager, and 36 were adults more than 16 years old. 
Patient age ranged from 0.2 to 90 years old (mean 29.6 
years; median 20 years), with 33 males and 32 females 
(Table 2). In 17 (26.2%) of all patients an additional virus 
was co-detected. Each of the genotypes showed co¬ 
detection except genotype C. The most frequent co¬ 
detected viruses were RSV and HRV. The age distributions 
in different genotypes differed significantly (One-way anal¬ 
ysis of variance, P = 0.0094). Genotype D was detected in 
patients with a broad age range (0.2—90 year old), although 
the majority (35 out of 51 cases) occurred in children and 
adults less than 50 years old. Genotype B was detected in 
one young adult (21 years old) with URTI in 2006, and in 
five children with LRTIs after 2010. Genotype C was de¬ 
tected in three older adults aged 58, 72 and 88 years with 
URTIs, whereas genotype E was only detected in five chil¬ 
dren less than 3 years of age (0.8—2.7 years old) with LRTIs. 

Discussion 

As an important human respiratory virus, the epidemic 
features of HCoV-OC43 at molecular level have not been 
well addressed. In this study, we describe the molecular 
epidemiological features of HCoV-OC43 in detail based on 
65 cases. Our results showed marked variations of HCoV- 
OC43 genotype prevalence from year to year, similar to 
that observed in other HCoVs. In line with previous 
reports, 1 genotypes B and C were detected before 
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Figure 2 Phylogenetic analysis of multiple gene regions of HCov-OC43 strains. A total of 23 gene regions are analyzed, including 
nspl to nspl 6, ns2a, HE, S, NS5ot, E, M and N of ten genomes of HCoV-OC43. The neighbor-joining method (Kimura’s two-parameter) 
was used to construct the trees with 1000 bootstrap values. 
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Figure 3 Recombination analyses of HCoV-OC43 genomes. Bootscan plot analysis of the whole genome sequences identified in 
this study in comparison to reference strains in genotypes B, C and D. Graphs were generated using 1783A/10, 2058A/10, 
3074A/12 and 3194A/12 as query sequences. The bootstrap value is 1000 for a window of 200 bp. 
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Figure 4 Yearly distribution of HCoV-OC43 genotypes during 2005—2012. The left and right vertical axis show the number of 
HCoV-OC43 cases and the detection rate among the recruited cases. 


2006, but disappeared after 2006. However, genotype B re- 
emerged in 2010 in our study, which has not been reported 
before. Genotype D, generated from recombination and 
first identified in 2004, was the dominant genotype for ep¬ 
idemics starting in 2007. This findings overlap with those of 
Lau et al., who reported that eight HCoV-OC43 strains de¬ 
tected between 2008 and 2011 all belonged to genotype 
D. 12 In our study, genotype D was not detected in 2011 
but in 2012, albeit at lower numbers (four out of seven 
HCoV-OC43 positive samples). It seems that immunity 


developed in the human population after the wide-spread 
of genotype D had blocked its epidemic as the overall prev¬ 
alence of genotype D showed decreased over time. Addi¬ 
tional analysis of the evolution of antigenic genes, 
particularly the S gene will help to further our understand¬ 
ing of the adaption of viral genotypes. 

Recombination is a common phenomenon among corona- 
viruses. A special random template switching mechanisms 
can be used during RNA replication. ,24 The high frequency 
of homologous recombination together with the high 


Table 2 Characteristics of HCoV-OC43 positive cases. 




Parameters 

Genotype 




Total 


B 

C 

D 

E 


Positive cases 

6 (9.2) a 

3 (4.6) 

51 (78.5) 

5 (7.7) 

65 

Age, years 

Range' 

0.6-21.0 

58.0-88.0 

0.2-90.0 

0.8-2.7 

0.2-90.0 

Mean/Median 

5.3/1.8 

72.7/72.0 

32.5/23.0 

2.5/1.1 

29.6/20.0 

Gender (M/F) 

4/2 

1/2 

26/25 

2/3 

33/32 

Diagnosis 

URTIs 

1 (16.7) 

3 (100) 

33 (64.7) 

0 

37 (56.9) 

LRTIs 

5 (83.3) 

0 

18 (35.3) 

5 (100) 

28 (43.1) 

Co-detected cases 

3 (50) 

0 

12 (23.5) 

2 (40.0) 

17 (26.2) 

1 virus 

3 (50) 

0 

8 (15.7) 

1 (20.0) 

12 (18.5) 

2 virus 

0 

0 

2 (3.9) 

1 (20.0) 

3 (4.6) 

3 virus 

0 

0 

2 (3.9) 

0 

2(3.1) 

Frequency of co-detected respiratory viruses 
RSVA/B 1 (16.7) 

0 

5 (9.8) 

2 (40.0) 

8 (12.3) 

HRV 

2 (33.3) 

0 

5 (9.8) 

0 

7 (10.8) 

IFVA 

0 

0 

2 (3.9) 

1 (20.0) 

3 (4.6) 

HPIVs 

0 

0 

3 (5.9) 

0 

3 (4.6) 

HBoV 

0 

0 

2 (3.9) 

0 

2 (3.1) 

hMPV 

0 

0 

1 (2.0) 

0 

1 (1.5) 


M, male. F, female. URTIs, upper respiratory tract infections. LRTIs, lower respiratory tract infections. RSV, respiratory syncytial virus. 
HRV, rhinovirus. IFVA, influenza viurs type A. HPIVs, human parainfluenza virus. hMPV, human metapneumovirus. HBoV, human 
bocaviurs. 

a Numbers in parentheses indicate the percentages of positive detections reported to the total number of positive samples. 
b P = 0.0094. 
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mutation rates of the genome may lead to the adaptation of 
CoVs and allow the generation of new strains and gen¬ 
otypes. For example, recombination has been reported 
to generate new genotypes and to contribute to the genetic 
diversity in HCoV-HKUl and -NL63, with recombination sites 
on nsp6/nsp7, nsp16/HE, and nsp3, and the S genes, respec¬ 
tively. Our work, which included a larger number of sam¬ 
ples over a longer surveillance period than previous studies, 
shows that a novel genotype E emerged in 2010. This high¬ 
lights again the role of recombination in the evolution of 
HCoV-OC43. Based on nucleotide identity comparison, phylo¬ 
genetic analysis of different genes, and boot scanning anal¬ 
ysis, genotype E might be generated from a recombination 
between genotypes B, C and D. Potential recombination sites 
may be at the junctions of nsp2/nsp3, nsp6/nsp7, nsp9/ 
nspIO, nsp12/nsp13, ns5oc/E and M/N gene. However, these 
observations need to be clarified based on more whole 
genome sequences of OC43. In addition, our results together 
with those of previous reports on the recombination analysis 
of HCoVs, indicate that the amplification of genes including 
at least nsp2/nsp3, nsp12/nsp13 (corresponding to pol gene) 
and the S and N genes is needed for genotyping and recom¬ 
bination analysis. 12 ’ 20 - 26 

The association of HCoV-OC43 genotypes with disease 
severity has not been well defined. A previous study found 
that among eight genotype D positive patients, seven were 
diagnosed with pneumonia. However, in our study, geno¬ 
type D showed no associations with severe symptoms, as 
most of the patients suffered from URTIs. This difference 
in results may be attributed to the studied cohort and num¬ 
ber of positive cases. However, host immune pressure in 
response to genotype D during a long epidemic period 
may also affect virulence. 

All cases of the novel genotype E and those of genotype 
B identified after 2010 were found in children younger than 
three years with LRTIs, but not detected in adults with 
LRTIs or URTIs. However, as the number of positive cases 
was limited, it is unclear whether the association of 
genotypes with LRTIs and special age groups is significant. 
This association may require further investigations for a 
larger number of samples. In addition, it should be further 
investigated whether the genetic configuration of genotype 
E allow it to spread rapidly, leading to the replacement of 
other genotypes such as genotype D. 

In summary, our results on the evolving genotypes of 
HCoV-OC43 and the emergence of a novel genotype E 
indicate that genotype shift may be one of the major 
ways for HCoV-OC43 to maintain its epidemic. Our findings 
provide insight into the evolution of HCoVs and its epide- 
micity, and can help inform CoV surveillance and control in 
humans and animals. 
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